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PRELIMINARY  INVESTIGATION  INTO  THE  EFFECT  OF 
I IIGII  TEMPERATURE  THERMAL  ANNEALS  Oil  THE  FIGURE  OF 
MERIT  OF  SILICON  GERMANIUM  ALLOYS 

ABSTRACT 

In  this  report  is  embodied  the  results  of  a  preliminary  theoretical 
and  experimental  investigation  into  the  effect  of  high  temperature 
thermal  anneals  on  the  thermoelectric  figure  of  merit  of  silicon 
germanium-gallium  phosphide  material. 

A  realistic  theoretical  model  for  the  silicon  germanium  alloy 
system  is  employed  to  investigate  the  effect  of  an  increase  in  the 
doping  level  on  the  electrical  power  factor  a2o  where  a  is  the  Seebeck 
roefficient  and  a  the  electrical  conductivity.  The  assumption  is  made 
that  the  addition  of  gallium  phosphide  to  silicon  germanium  alloy 
serves  only  to  increase  the  dopant  solubility  and  hence  the  carrier 
concentration.  The  theoretical  analysis  indicates  that  the  reported 
factor  of  1.6  increase  in  the  carrier  concentration  of  heat  treated 
silicon  germanium-gallium  phosphide  material  compared  to  conventional 
silicon  germanium  alloy  would  result  in  an  increase  in  the  electrical 
power  factor  of  between  12-14%  at  room  temperature  rising  to  18-22% 
at  1000K. 

Comparison  of  measured  transport  data  on  silicon 
germanium-gallium  phosphide  material  and  a  zone  levelled  silicon 
germanium  alloy  confirms  that  substantial  improvements  in  the 
electrical  power  factor  can  be  obtained  by  subjecting  the  silicon 
germanium-gallium  phosphide  material  to  high  temperature  thermal 
anneals,  while  the  total  thermal  conductivity  remains  relatively 
unchanged . 

Although  there  are  differences  in  the  increases  in  the  power 
factor  obtained  by  different  researchers.  The  results  of  this 
programme  of  work  support  the  conclusion  that  the  electrical  power 
factor  and  hence  the  figure  of  merit  of  silicon  germanium-gallium 
phosphide  can  be  significantly  improved  by  high  temperature  thermal 
anneals. 
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Ill  General  Introduction 

Silicon  germanium  alloys  are  established  semiconductor  material 

for  high  temperature  applications  and  they  have  been  successfully 

employed  in  the  fabrication  of  thermocouples  of  radioisotopic  powered 

generators  which  provided  onboard  electrical  power  to  the  Lincoln 

Experimental  Satellites  LES  8  and  LES9  and  the  voyager  spacecrafts1. 

Materials  based  upon  the  silicon  germanium  alloys  have  also  been1 

chosen  for  use  in  NASA's  SP-100  project  in  which  orbiting  nuclear 

reactor  powered  thermoelectric  generators  will  produce  up  to  100 

2 

kilowatt  of  electrical  power  . 

In  space  applications  weight  is  a  major  consideration  and  over 

the  past  ten  years  or  so  effort  have  been  made  to  reduce  the  weight 

of  the  thermoelectric  conversion  system  by  improving  the  performance 

of  the  thermocouple  material.  The  material's  performance  is 

conveniently  expressed  in  terms  of  its  figure  of  merit  Z  =  z2a/X  where 

a  is  the  Seebeck  coefficient,  a  the  electrical  conductivity  and  X  the 

thermal  conductivity.  The  thermal  conductivity  consists  essentially  of 

two  parts,  a  lattice  (phonon)  contribution  A^  and  an  electron  (hole) 

contribution  X^.  The  lattice  thermal  conductivity  is  the  major 

component,  even  in  very  heaviiy  doped  thermoelectric  materials. 

Effort  has  focussed  on  improving  Z  and  hence  increasing  the 

conversion  efficiency  by  reducing  the  lattice  component  of  the  thermal 

conductivity.  Substantial  reductions  in  A.  have  been  achieved 

L  3 

through  the  use  of  very  small  grain  size  material  and  through  the 
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introduction  of  additives  such  as  gallium  phosphide  .  However,  all 
published  information  indicates  that  although  a  reduction  in  A^  can  be 
obtained  any  improvement  in  Z  is  minimal  due  to  an  accompanying 
degradation  in  electrical  properties  -  primarily  a  reduction  in  carrier 
mobility . 

Recently  attention  has  turned  to  improving  the  electrical  power 
factor  (EPE).  It  has  been  reported  that  the  EPF  of  silicon  germanium 
-  gallium  phosphide  can  be  significantly  increased  by  subjecting  the 
material  to  high  temperature  anneals.  This  increase  in  the  EPF  is  not 
iccompanied  by  a  significant  chance  in  the  thermal  conductivity  and 
results  ir:  an  improvement  in  the  thermoelectric  figure  of  merit'5.  A 
suggested  mechanism  responsible  for  the  reported  increase  in 
performance  is  that  the  gallium  phosphide  enhances  the  dopant 
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solubility  hence  permitting  the  carrier  concentration  to  approach 
closer  to  its  optimum  valued  These  developments  lead  to  the  present 
programme  of  research. 

IV  Objectives 


The  programme  of  work  proceeded  in  two  stages,  the  first  one 
theoretical  and  the  second  one  experimental.  In  the  first  stage  the 
realistic  theoretical  model  for  silicon  germanium  alloy,  which  had  been 
developed  for  a  previous  U.S.  Army  contract7  would  be  used  to 
explore  the  potential  improvement  in  the  electrical  power  factor  and 
hence  thermoelectric  figure  of  merit  which  accompanies  an  increase  in 
dopant  solubility. 

The  experimental  programme  of  work  involves  annealing 
specimens  of  silicon  germanium-gallium  phosphide  at  high  temperatures 
and  measuring  the  Scebeck  Coefficient  and  electrical  conductivity. 
Finally  an  assessment  would  be  made  of  the  effect  of  high  temperature 
anneals  on  the  electrical  power  factor  of  silicon  germanium-gallium 
phosphide  ana  possible  mechanisms  suggested  which  are  responsible 
for  any  observed  increase  in  the  materials  performance. 

V  Theoretical  Programme  of  Work 


1 .  Introduction 


Silicon  germanium  alloys  cannot  be  described  as  narrow 

band  semiconductors,  however  the  high  level  of  doping  normally 

employed  in  thermoelectric  materials  necessitates  the  inclusion  of  a 

deviation  from  the  usually  assumed  parabolic  bands.  The 

non-parabolicity  of  the  bands  is  interpreted  in  terms  of  a  parameter 

6  =  knT/E  where  kn  is  the  Baltzmann  constant  and  T  is  the  absolute 
ti  g  d 

temperature.  E  t  is  the  energy  band  gap.  Non-parabolicity  has  a 
9 

strong  influence  on  the  electrical  properties  at  high  carrier 
concentration.  Moreover,  these  effects  are  quantitatively  dependent 
on  the  scattering  mechanicam.  The  development  of  a  comprehensive 
and  realistic  theoretical  model  which  provides  a  description  of  the 
electrical  properties  of  the  silicon  germanium  alloy  system  is  difficult 
and  best  approached  in  stages,  commencing  with  a  working  model  and 
progressing  to  more  realistic  models  of  increased  complexity. 


At  room  temperature  and  a  little  above,  acoustic  phonon 
scattering  is  the  dominant  charge  carrier  scattering  mechanism .  In 
the  model  adopted  we  have  considered  a  multivalued,  non  parabolic: 
band  structure  but  excluded  the  effect  of  intervally  scattering, 
intervctlly  scattering  requires  the  presence  of  high  energy  phonons 
which  are  not  excited  at  around  room  temperature 
2.  Theoretical  model 


It  is  convenient  to  work  in  terms  of  reduced  electrical 
conductivity  a1  and  reduced  Seebeck  coefficient,  a1  where 
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with  e j  the  accoustic  deformation  potential,  C  ..  the  longitudinal 
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elastic  constant  Ny  the  number  of  equivalent  energy  valleys,  mc  the 
conductivity  (intertial  effective  mass)  and  A  ^  the  lattice  thermal 
conductivity.  The  generalised  Fermi  integrals  nL^  which  appears  in 
the  electronic  transport  coefficient  expressions  are  given  by: 


lr  m 


U) 


-  ff 
a 


n 


]r?n  [n(l  +  0n)]m  (1  +  20n) 


dn 


f  is  the  Fermi  distribution  function,  n  =  E/kT,  8  =  knT/E 

b  g 

The  carrier  concentration  n  is  related  to  the  Fermi  level 

c 

nc  =  Nv  (2md*s  kBT)V2  oLy 
3n2 


by: 


The  density  of  states  yo(Ej  of  a  single;  valley  is  given  by: 


r  /-HI Af  \  i  /  / 

gs(E)  =  m  a-’  ft  -  di 

h2 


where  E  is  the  carrier  energy  measured  with  respect  to  t Mo  band 

* 

edge  and  is  the  density  of  states  effective  mass  for  a  single 

valley. 

If  there  are  N  eauivalent  valleys,  the  total  density  of 
states  g  (E)  is  given  by: 

gt(E)  =  Nv(gs(E) 

This  can  then  be  incorporated  in  our  model  by  defining  a  total 

* 

density  of  states  effective  mass  ni^  such  that 


gt(E)  -  4 77  f  1 
h2  J 


3'  2 


E*dE 


This  gives  m^jt  =  Ny2/3m^s  =  N  v  2/  3(m  ,  m?  rn-j) v  3 

where  ntj,  m2  and  arc  the  components  of  the  mass  tenser  along  the 
principal  axis. 

The  effective  mass  values  for  the  electronic  states  in  the  alloys 
are  about  the  same  as  for  corresponding  states  in  pure  crystals, 
consequently  initially  a  simple  linear  compositional  dependence  of  the 
effective  mass  was  employed  using  effective  mass  values  given  in 
Table  1,  making  no  distinction  between  the  density  of  states  and 
conductivity  effective  masses. 


Ge 

Si 

m.  / 1 
t  0 

0.08 

0.19 

m.  /m 

L  0 

1.58 

0.92 

m  ,/m 
d  0 

0.216 

0.32 

m  /  m 

0.117 

0.26 

Table  1  -  Effective  mass  values 
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The  values  for  and  E^  are  taken  to  be  1.7  x  10  Nm  and 

6.2eV  respectively.  The  calculations  have  been  carried  out  for  70% 
and  90%  silicon  in  view  of  the  expected  usefulness  of  silicon  rich 


alleys  in  ''ctua1  device  applications.  A  slight  improvement  in  the 
model  can  be  obtained  by  considering  the  density  of  states  and 
conductivoly  effective  masses  separately  for  silicon  and  germanium  and 
obtaining  the  corresponding  values  for  the  alloy  composition  by  linear 
interpolations. 

3 .  Theoretical  Resui ts 

In  the  theoretical  analysis  it  is  assumea  that  the  normal 
material  doping  level  is  well  below  the  maximum  required  to  optimise 
the  figure  of  merit  and  that  no  significant  changes  in  the  effective 
mass  on  the  other  bond  structure  parameters  of  silicon  germanium 
alloy  accompany  the  addition  of  gallium  phosphide.  The  theoretical  Z 
versus  carrier  concentration  curve  remains  unaltered  as  shown 
schematically  in  Figure  1.  With  the  carrier  concentration  n^  being 
determined  by  the  solubility  limit  of  the  dopant.  The  addition  of 
gallium  phosphide  enables  the  carrier  concentration  to  increase  to  a 
value  n^  which  approaches  closer  to  the  maximum  of  the  Z  versus 
carrier  concentration  curve.  The  theoretical  model  described  in 
section  V.  2.  has  been  employed  in  obtaining  the  dependence  of  the 
electrical  power  factor  on  carrier  concentration  for  different  alloy 
compositions  at  room  temperature  (Figure  2)  and  at  lOOOK  (Figure  3). 
Using  Figures  2  and  3  it  is  possible  to  estimate  the  increase  in  the 
electrical  power  factor  which  corresponds  to  the  factor  of  1 . 6  increase 
in  carrier  concentration  measured  at  500°C,  which  accompanies  high 
icmpc-ature  anneals  . 

The  results  are  presented  in  Tables  2  and  3. 

Composition  (%bi)  Carrier  Concentration  E.P.F.  %  increase  in 


,  -3. 

(cm  ) 

(W/mK2) 
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70 

io19 

40x1C~3 

1 ,6x1019 

45x1 0~3 

12.5 

90 

io19 

31  xK 

1  9 

1.6x10 

40x1 03 

14.3 

Table  2. 

Theoretical  increase  in 

the  electrical  power 

factor  at 

room  temperature. 
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2  2 

I  cb!t>  3.  Theoretical  increase  in  the  electrical  power  factor  at 
1  000  K  . 

4.  Piscussion 

The  theoretical  increase  in  the  electrical  power  factor  which 
accompanies  an  increase  in  the  carrier  concentration  by  a  factor  of 
i.6  range  from  12-14;  at  room  temperature  rising  to  18-22%  at  1000K. 
At  500  C  this  increase  is  about  half  the  reported  value  of  42%.  This 
simple  model  which  assumes  that  the  addition  of  CaP  only  causes  an 
increase  in  the  dopant  solubility  is  unable  to  account  for  the  reported 
increase  in  the  electrical  power  factor.  Either  the  measured  increase 
is  incorrect  or  the  theoretical  model  requires  refining  to  take  into 
account  changes  in  band  structure  and  density  of  states  which  may 
accompany  the  addition  of  GaP  to  silicon  germanium  alloy. 

VI.  Experimental  Work 

1.  Introduction 

The  ther  rnoelectr  ic  properties  of  heavily  doped  silicon 
germanium  alloys  are  notoriously  dependent  upon  the  history  of  the 
material  and  the  same  is  probably  true  for  silicon  germanium-gallium 
phosphide.  Considerable  differences  in  the  magnitude  in  the  increase 
in  the  electrical  power  factor,  which  accompanies  high  temperature 

5 

anneals,  have  been  reported  by  different  laboratories  .  Evidently 
additional  experimental  data  on  these  materials  is  required  and  this 
prompted  the  following  limited  programme  of  experimental  work. 


Procedure 


Thermoelectric  transport  property  measurements  wore  made 
ut".  a  sample  of  silicon  germanium-pallium  phosphide  material  obtained 
from  the  Jet  Propulsion  Labora  tor  ies ,  Pasadena.  An  ingot  shaped 
specimen  2mm  x  2mm  x  6mm  was  cut  from  a  pressed  pellet.  Its 
density  measured  using  Archemedes  principle  was  2.95  umcm  ^ . 
X-ray  analysis  indicated  that  the  alloy  composition  was  about 
SigE-Ge^t-.  The  specimen  was  head  treated  in  air  at  1  500 K  for  two 
hours  and  air  quenched.  It  was  subsequently  heated  in  air  at  lower 
temperatures  in  the  range  1  400K  to  HOOK.  The  heat  treatment  was 
periodically  interrupted  by  quenching  tne  specimen  in  air  to  room 
temperature.  An  oxide  layer  which  formed  on  the  specimen  only  after 
the  initial  high  temperature  anneal  was  removed  with  HF(403). 

Prior  to  measurements  the  surface  of  the  specimen  was  polished  and 
washed  in  distrilled  water.  Electrical  resistivity  and  Seebeck 
coefficient  measurements  were  made  from  room  temperature  to  about 
330K  . 


3 .  Results 

The  measured  Seebeck  coefficient  and  electrical  resistivity 
for  a  silicon  germanium-galliurn  phosphide  specimen  are  presented  in 
Figure  4,  5  and  6.  In  Figure  7  are  displayed  the  electrical  power 
factors  as  measured  at  different  laboratories  together  with 
measurements  made  some  time  agao  at  Cardiff  on  a  zone  levelled 
specimen  of  silicon-germanium  alloy  obtained  from  the  David  Sarnoff 
Centre  (RCA). 

4 .  Discussion 

The  silicon  germanium-gallium  phosphide  alloy  which  had 
been  subjected  to  a  high  temperature  thermal  anneal  exhibited  an 
increased  electrical  figure  of  merit  compared  to  standard  silicon 
germanium  materials.  At  the  highest  temperature  of  measurement  i.e. 
about  830K  the  increase  in  the  electrical  power  factor  was  .  This 
compared  favourably  with  the  results  reported  by  TECO  but  was 
considerably  less  than  the  increase  reported  by  JPL. 


VII  General  Discussion  and  Conclusions 

The  experimental  results  obtained  support  the  proposition  that 
the  electrical  power  factor  can  be  significantly  increased  by 
subjecting  silicon  germanium-gallium  phosphide  material  to  high 
temperature  thermal  anneals.  The  measurements,  although  only 
undertaken  on  one  sample  indicates  that  the  percentage  increase  in 
electrical  power  factor  is  somewhat  less  than  that  obtained  at  the  Jet 
Propulsion  Laboratory  and  about  the  same  order  as  reported  by 
TECO.  A  theoretical  analysis  based  upon  a  realistic  model  for  silicon 
germanium  in  which  it  is  assumed  that  the  presence  of  gallium 
phosphide  serves  only  to  increase  the  solubility  of  the  dopant,  would 
account  for  this  level  of  improvement.  However,  the  model  adopted  is 
unable  to  explain  the  much  greater  improvement  reported  by  workers 
at  JPL. 
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FIGURE  OF  MERIT  Z  — ~  FIGURE  OF  MERIT  Z 


FIG. I  (a)  FIGURE  OF  MERIT  versus  n  (schematic). 
SIMPLE  MODEL:  NO  BAND  STRUCTURE  CHANGES  WITH 
GaP  DOPING  :  ONLY  SOLUBILITY  INCREASE. 


FIG.  1(b)  FIGURE  OF  MERIT  versus  n  (schematic). 
IMPROVED  MODELiBAND  STRUCTURE  CHANGES 
ACCOMPANY  SOLUBILITY  INCREASE. 


{wo  uh{oid )  AJ.IAIJLSIS3H  1VDIHA3H13 


temperature  (K) 

FIG.5  ELECTRICAL  RESISTIVITY  versus  TEMPERATURE 


IMPROVED  JPL  SAMPLE 


FIG  7  COMPARISON  OF  THE  ELECTRICAL  POWER  FACTOR 
DATA  ON  SiGe-GaPAS  A  FUNCTION  OF  TEMPERATURE 


